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Abstract 

This paper presents necessary and sufficient optimality conditions for problems of the 
fractional calculus of variations with a Lagrangian depending on the free end-points. The 
^ ' fractional derivatives are defined in the sense of Caputo. 
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1 Introduction 



Fractional calculus is one of the generalizations of the classical calculus and it has been 
used successfully in various fields of science and engineering — see, e.g., [IT ^fT^ fTo ^ fTU lEZOj 
[231 12H ISSl EH HEB EH ES]. In recent years, there has been a growing interest in the area of 
fractional variational calculus and its applications plMn ^ [T ^[T6iri8l [22 l l28 p 2] . Applications 
include classical and quantum mechanics, field theory, and optimal control. In the papers 
cited above, the problems have been formulated mostly in terms of two types of fractional 
derivatives, namely Riemann-Liouville and Caputo. The natural boundary conditions for 
fractional variational problems, in terms of the Riemann-Liouville and the Caputo derivative, 
are presented in [U[2]. Here we develop further the theory by proving necessary optimality 
conditions for more general problems of the fractional calculus of variations with a Lagrangian 
that may also depend on the unspecified end-points y(a), y(b). More precisely, the problem 
under our study consists to minimize a functional which is defined in terms of the Caputo 
derivative and having no constraint on y(a) and/or y(b). The novelty is the dependence of 
the integrand L on the free end-points y(a), y(b). This class of problems is motivated by 
applications in the field of economics [10J. 

The paper is organized as follows. Section[2]presents the necessary definitions and concepts 
of the fractional calculus; our results are formulated, proved, and illustrated through an 
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example in Section [3l Main results of the paper include necessary optimality conditions with 
the new natural boundary conditions (Theorem 13. ip that become sufficient under appropriate 
convexity assumptions (Theorem I3.3j) . 



2 Fractional Calculus 

In this section we review the necessary definitions and facts from fractional calculus. For 
more on the subject we refer the reader to [2T | [2l ^ [3"T | l3"3] . 

Let / G Li([a, b\) and < a < 1. We begin with the left and the right Riemann-Liouville 
Fractional Integrals (RLFI) of order a of function / which are defined as: the left RLFI 

alZfix) = (\x - tr^f^dt, x e [a, b], (1) 

the right RLFI 

xlbf(^) = TT7~T / (t~ x) a - l f(t)dt, x G [a, b], (2) 



r(a) Jx 

where T(-) represents the Gamma function. Moreover, a I®f = x Ij!f = / if / is a continuous 
function. 

Let / G AC([a,b\), where AC([a,b]) represents the space of absolutely continuous func- 
tions on [a, b\. Then using equations ([I]) and ([2]), the left and the right Riemann-Liouville and 
Caputo derivatives are defined as: the left Riemann-Liouville Fractional Derivative (RLFD) 

a D a J{x) = Y^Taj-^ £( x ~ trV(t)dt = ^Jl- a f{x), x G [a, b], (3) 
the right RLFD 

the left Caputo Fractional Derivative (CFD) 

aD«f(x) = f7 ^y J\x - ty a j t f(t)dt = Jl-±f{ x ), x g [a, b], (5) 
the right CFD 

xDtf{x) = Tirh) iy ~ x) ~ a jt fit)dt = xI ^ a x G [a ' 6] > (6) 

where a is the order of the derivative. 

The operators {I])-© are obviously linear. We now present the rule of fractional integra- 
tion by parts for RLFI (see for instance [S]). Let < a < 1, p > 1, q > 1, and 1/p+l/q < l+at. 
If g G L p ([a,b]) and / G L q ([a,b]), then 

b g(x) a I«f(x)dx = I" f(x) x I?g(x)dx. (7) 
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In the discussion to follow, we will also need the following formulae for fractional integrations 
by parts: 

* g{x) c a D«f{x)dx= f(x) x I 1 b - a g(x)\ X x Z b a + f f(x) x D^g(x)dx, 
/ g(x)%D?f(x)dx= -f(x) a I 1 - a g(x)\ X x Z a + / f(x) a D a x g{x)dx. 

J a J a 

They can be derived using equations ©-flSJ), the identity ((7J) and performing integration 
by parts. 

3 Main Results 

Let us consider the following problem: 

rb 

J(y) = J L(x,y(x),^D^y{x),^D^y(x),y(a),y(b))dx — > extr 

= y a ), (y(b) = y b ). 

Using parentheses around the end-point conditions means that the conditions may or may 
not be present. We assume that: 

(i) L( v , .,.,.,.) eC'fKtlxl^); 

(ii) x — > d^L(x, y(x), < ^D"y(x), x D b y(x),y(a),y(b)) has continuous right RLFI of order 1—a 
and right RLFD of order a, where a € (0, 1); 

(hi) x — > diL(x, y(x),a D"y(x),^ D®y(x),y(a),y(b)) has continuous left RLFI of order 1 — /3 
and left RLFD of order f3, where /3 € (0, 1). 

Remark 1. We are assuming that the admissible functions y are such that ^D x y(x) and 
x D b y(x) exist on the closed interval [a,bj. 

Along the work we denote by d{L, i = 1, . . . , 6, the partial derivative of function L(-, •, ■, •) 
with respect to its ith argument. 

3.1 Necessary Optimality Conditions 

Next theorem gives necessary optimality conditions for the problem ([9]). 

Theorem 3.1. Let y be a local extremizer to problem ©. Then, y satisfies the fractional 
Euler- Lagrange equation 

d 2 L(x, y(x)^D^y(x)^D b y(x),y(a),y(b))+ x D^d 3 L(x, y{x)° D«y{x)° D^y(x),y(a),y(b)) 

+ a D%L(x,y(x)CD«y(x), x 1 D b 3 y(x),y(a),y(b)) = (10) 
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for all x G [a,b]. Moreover, if y(a) is not specified, then 
/ d 5 L(x, y(x),^D^y(x),^D^y(x),y(a),y(b)) dx 

J a 

= ' x ll- a d 3 L(x, y{x)° D«y(x)C D Py(x),y(a),y(b)) 

- a I x ^d 4 L(x,y(x)^D^y(x)^D^ b y(x),y(a),y(b)) 

if y(b) is not specified, then 

/ d G L(x, v(x),2n%y(x),%D%y(x),y(a),y(b)) dx 

J a 

a I^d 4 L(x, y{x), c a D«y{x), c x D P b y{x),y{a),y(b)) 

- x ll- a d 3 L(x, y{x)° D%y(x)CDPy(x),y(a),y(b)) 



(11) 



x=b 



(12) 



Proof. Suppose that y is an extremizer of J . We can proceed as Lagrange did, by considering 
the value of J at a nearby function y = y + eh, where e £ M is a small parameter, h is an 
arbitrary admissible function. We do not require h(a) = or h(b) = in case y(a) or y(b), 
respectively, is free (it is possible that both are free). Let 



<p(e) = I L(x,y(x)+eh(x),^D^(y(x)+£h(x)),^D^(y(x)+eh(x)),y(a)+eh(a),y(b)+eh(b))dx 

J a 

Since and x D b are linear operators, it follows that 

c a D«(y(x) + eh{x)) = ^D^y(x) + e^D^h(x) 
C x D P b {y(x) + eh(x)) = C x D P b y{x) + e C x D P b h(x). 



A necessary condition for y to be an extremizer is given by 
= 



d(p 
de 



e=0 



d 2 L(- ■ ■ )h(x)+d 3 L(- ■ ■ ) c a D«h{x)+d±L{- ■ ■ )^h(x)+d 5 L(- ■ ■ )h(a)+d 6 L(- ■ ■ )h(b) 



dx 



= 0, (13) 

where (• • • ) = (x,y(x),^ D x y(x), x D b y(x),y(a),y(b) S j . Using formulae (jSJ) for integration by 
parts, the second and the third integral can be written as 



d 3 L(.-.f a D«h(x)dx= I h(x) x D^d 3 L(- ■ -)dx + x I^ a d 3 L(. ■ .)h(x)\lZ u a , 



i x=b 



rb pb 

I d 4 L(- ■ ■ )^D p b h{x)dx = I h{x) a D^diL(- --)dx- a I l -^d 4 L{- 

J a J a 



)h(x) 



x=b 



(14) 
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Substituting (JTU) into (|T3j). we get 

rb 



/ [5 2 L(- • • ) + x D%d 3 L(- ■■ )+ a D%L(- ■ • )1 fc( 5 

J a 



+ xl 



1-a 



d 3 L(---)h(x)\ X x Z b a -aI x ^d 4 L(...)h(x) 

b 



x=b 



+ 



(d 5 L(- ■ ■ )h(a) + d 6 L(- ■ ■ )h(b)) dx = 0. (15) 



We first consider functions h(t) such that h{a) = h(b) = 0. Then, by the fundamental lemma 
of the calculus of variations, we deduce that 



d 2 L(- ■■)+ x D£d 3 L(- • • ) + a D%L(- • • ) = 



(16) 



for all x € [a, b\. Therefore, in order for y to be an extremizer to the problem @, y must be 
a solution of the fractional Euler-Lagrange equation. But if y is a solution of f 1 1 6 [) . the first 
integral in expression (|15p vanishes, and then the condition (|13p takes the form 



Kb) 



d 6 L(- ..)dx- [j x - p d±L{- ■■)- x ll- a d 3 L(- ■ ■ 



\x=b 



+ h(a) 



d$L(- ■ ■ )dx 



0. 



If y(a) = y a and y(b) = yi, are given in the formulation of problem ([9]), then the latter equation 
is trivially satisfied since h{a) = h(b) = 0. When y(a) is free, then 



J a 



d 5 L(- ..)dx- [xil- a d 3 L{- ■■)- a i x -?d 4 L(- ■ ■ ) 

when y(b) is free, then 

d 6 L(- --)dx- \ a I^d A L(- • • ) - xll~ a d 3 L{- ■ ■ ) 



0, 



since h(a) or h(b) is, respectively, arbitrary. 



□ 



Remark 2. Conditions (1101 ) -( 1121) are only necessary for an extremum. The question of 
sufficient conditions for an extremum is considered in Subsection \3.£\ 



In the case L does not depend on y(a) and y(b), by Theorem 13. II we obtain the following 
result. 

Corollary 1 (Theorem 1 of [2]). If y is a local extremizer to problem 

,-b 



J(y) 



L {x,y(x),^D x y(x), x D^y(x))dx — > extr. 



then y satisfies the fractional Euler-Lagrange equation 
d 2 L(x, y(x), C a D«y{x), %D P h y{x)) + x D^d 3 L(x, y(x), C a D«y{x), C x D P b y{x)) 



+ a D%L(x, y(x), c a D a x y(x), c x D p h y{x)) = 



5 



for all x G [a,b]. Moreover, if y(a) is not specified, then 

4- a d 3 L(x, y{x), c a D«y{x)CD P b y(x)) - a I x ^d 4 L(x, y{x)°L%y{x), °D P h y{x)) 
if y(b) is not specified, then 

' a I*- p d4L(x,y{x)°I%y(x)°Dty(x)) -^^{x, y(x), C a D«y{x), C X D? y{x)) 



0, 



= 0. 

x=b 



We note that the generalized Euler-Lagrange equation contains both the left and the right 
fractional derivative. The generalized natural conditions contain also the left and the right 
fractional integral. Although the functional has been written only in terms of the CFDs, neces- 
sary conditions (I10p -( fl2l) contain Caputo fractional derivatives, Riemann-Liouville fractional 
derivatives and Riemann-Liouville fractional integrals. 

Observe that if a goes to 1, then the operators ^D x and a D x can be replaced with 4^ and 
the operators X D% and x D b can be replaced with — ^ (see [31]). Thus, if the xD^y term is 
not present in Q, then for a — > 1 we obtain a corresponding result in the classical context 
of the calculus of variations [10] (see also [25^ Corollary 1]). 

Corollary 2. If y is a local extremizer for 

J(y)= / L (x,y(x),y r (x),y(a),y(b))dx — > extr 



{y(a) = Ua), (y( b ) = Vb), 

then 

—d 3 L(x, y(x),y'(x), y(a),y(b)) = d 2 L(x, y(x), y'(x),y(a),y(b)) 
for all x G [a,b]. Moreover, if y(a) is free, then 

d 3 L(a,y(a),y'(a),y(a),y(b)) = / d s L(x,y(x),y'(x),y(a),y(b))dx) 

J a 

if y(b) is free, then 

d 3 L(b,y(b),y'(b),y{a),y(b)) = - d 6 L(x,y(x),y\x),y(a),y(b))dx. 

J a 

3.2 Sufficient Conditions 

In this section we prove sufficient conditions that ensure the existence of minimum (maxi- 
mum). Similarly to what happens in the classical calculus of variations, some conditions of 
convexity (concavity) are in order. 

Definition 3.2. Given a function L, we say that L(x,y, z,t,u,v) is jointly convex (con- 
cave) in (y, z,t,u,v), if d{L , i = 2,...,6 ; exist and are continuous and verify the following 
condition: 

L(x, y + y\, z + z\, t + 1 1} u + ui, v + v{) - L(x, y, z, t, u, v) 

> {<)d 2 L(*)yx + d 3 L(.) Zl + d 4 LC)ti + d 5 L(.) Ul + d e L(m) Vl 

for all (x,y, z,t,u,v), (x,y + y±,z + Z\,t + t\,u + u±,v + v\) G [a,b] x M 5 ; where (•) = 
(x,y,z,t,u,v). 
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Theorem 3.3. Let L(x,y, z,t,u,v) be jointly convex (concave) in (y, z,t,u,v). If yo satisfies 
conditions (|10p ~ (|12p . then y$ is a global minimizer (maximizer) to problem ([9]). 

Proof. We shall give the proof for the convex case. Since L is jointly convex in (y, z, t, u, v) 
for any admissible function yo + h, we have 



J{yo + h)- J{y ) 

rb r 



L(x, y Q (x) + h(x)°D2(yo(x) + h{x))^D p b (y Q {x) + h(x)),y (a) + h(a),y (b) + h(b)) 



L(x, yo(x),aD^y (x),^D^y Q (x),yo(a),y (b)) 



dx 



> 



d 2 L(*)h(x) + d 3 L(*)C D°h(x) + d 4 L(*)^D£h(x) + d 5 L(*)h(a) + d 6 L(*)h(b) 



dx 



0. 



where (★) = (x, yo(x), % D%yo(x), x yo(x), yo(a), yo(b)j . We can now proceed analogously 
to the proof of Theorem 13.11 As the result we get 

J(yo + h)- J(y ) > / <9 2 L(*) + x D£d 3 L(*) + a DPd A L{*) h(x) 

J a 

+ h{b) y'd 6 L(*)dx - [ a J^9 4 iW - l^j 

+ h(a) y'd 5 L(*)dx - [jj'-^lW - „t^ 4 LW 
Since yo satisfy conditions (fTUj) ~ (fT2l . we obtain J~(yo + h) — J(yo) > 0. 
3.3 Example 

We shall provide an example in order to illustrate our main results. 
Example 1. Consider the following problem: 

J(V) = \J [{%D a x y(x)f + 7 y 2 (0) + A(y(l) - l) 2 ] dx — > min 



□ 



(17) 



where 7, A € R + . For this problem, the generalized Euler-Lagrange equation and the natural 
boundary conditions (see Theorem \3.1\) are given, respectively, as 



Df {°D«y{x)) = 0, 



[\y(0)dx = x l{- a {CD«y(x)) \ x=0 , 
Jo 

A(y(l) - l)dx = - x lt a ($D%y(x)) \ x=1 . 



(18) 
(19) 
(20) 



Note that it is difficult to solve the above fractional equations. For < a < 1 numerical 
method should be used. When a goes to 1 problem (|17p tends to 



J(y) 



1 



[(y'(x)) 2 + 1 y 2 (0) + X(y(l)-l) 2 ]dx 



mm 



(21) 
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and equations f|18[) — (j20|) could be replaced with 



y"(x) = 0, 



(22) 



72/(0) = y'(0) 
A(y(l) -!) = -» 




(23) 
(24) 



Solving equations (|22p - (|24p we obtain that 



y{x) = 



7A 



x + 



A 



7A + A + 7 



7A + A + 7 



is a candidate for minimizer. Observe that problem (|17p satisfies assumptions of Theorem \3.3[ 
Therefore y is a global minimizer to problem (|21h . 
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